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Abstract 

The General Antiparticle Spectrometer (GAPS) is a novel approach for indirect dark matter searches that exploits 
cosmic antiparticles, especially antideuterons. The GAPS antideuteron measurement utilizes distinctive detection 
methods using atomic X-rays and charged particles from the decay of exotic atoms as well as the timing and stopping 
range of the incoming particle, which together provide excellent antideuteron identification. Prior to the future balloon 
experiment, an accelerator test and a prototype flight were successfully conducted in 2005 and 2012 respectively, in 
order to verify the GAPS detection concept. This paper describes how the sensitivity of GAPS to antideuterons was 
estimated using a Monte Carlo simulation along with the atomic cascade model and the Intra-Nuclear Cascade model. 
The sensitivity for the GAPS antideuteron search obtained using this method is 2.0 xl0“^ [m“^s“^sr“^(GeV//t)“^] for 
the proposed long duration balloon program (LDB, 35 days x 3 flights), indicating that GAPS has a strong potential 
to probe a wide variety of dark matter annihilation and decay models through antideuteron measurements. GAPS is 
proposed to fly from Antarctica in the austral summer of 2019-2020. 
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1. Introduction 

1.1. Dark Matter and WlMPs 

The recent result by the Planck experiment m shows 
that 68% of our universe is composed of dark energy, 
27% is dark matter and 5% is baryonic matter. The na¬ 
ture and origin of dark energy and dark matter, however, 
are still unknown. This is one of the great cosmologi¬ 
cal problems of the 21st century. The existence of dark 
matter was postulated by Fritz Zwicky in 1933 from the 
observation of the Coma galaxy cluster O • Recent ob¬ 
servations, such as gravitational lensing in the Bullet 
Cluster (two colliding clusters of galaxies), also indi¬ 
cate the existence of dark matter 0. Weakly Interacting 
Massive Particles (WIMPs) are among the theoretically 
best-motivated candidates in the variety of dark mat¬ 
ter models. The lightest supersymmetric partner (LSP) 


* Current institution: SLAC National Accelerator Labora- 
tory/Kavli Institute for Particle Astrophysics and Cosmology, Menlo 
Park, CA 94025, USA 

Email address: tsuguoOslac. stanford.edu (T. Aramaki) 
Preprint submitted to Astroparticle Physics 


in supersymmetric (SUSY) theories such as neutrali- 
nos HID and right-handed sneutrinos mn) as well as 
Kaluza-Klein particles, such as right-handed neutrinos 
(Lzp) m in extra dimension theories are examples of 
popular WIMP candidates. Gravitinos, known as Super 
Weakly Interacting Massive Particles (SuperWIMPs), 
are also favored dark matter candidates in SUSY the¬ 
ories 121. 

1.2. Antideuterons as Dark Matter Search 

Dark matter searches are usually categorized into 
three types: particle colliders, direct searches, and in¬ 
direct searches. The direct and indirect searches will 
measure the relic WIMPs in our universe, while the par¬ 
ticle collider will try to create WIMPs. The detection 
methods and the background models for each search are 
different, but also complementary, helping to illuminate 
the nature of dark matter. 

Indirect dark matter searches measure dark matter 
co-annihilation and decay products, such as electrons, 
positrons, gamma-rays, antiprotons and antideuterons. 
The General Antiparticle Spectrometer (GAPS) is a 
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Figure 1: Antideuteron flux at the top of the atmosphere. The solid 
blue and purple lines are neutralino and gravitino LSPs with ~ 30 
GeV O and m^ii ~ 50 GeV (9l respectively, while the dashed green 
line is for LZP with ~ 40 GeV . The red solid line represents the 
secondary/tertiary flux due to the cosmic-ray interactions ll2lll3lfT^ . 
The expected GAPS antideuteron sensitivity (~ 99% CL) as well as 
the AMS-02 sensitivity (Ts] [Tb] [TT] [Ts) and the current upper limit 
on the antideuteron flux obtained by BESS EH are also shown (see 
Section [33) . 


novel approach for indirect dark matter searches that ex¬ 
ploits cosmic antiprotons and antideuterons (TOl [TTl . 

Antideuteron production in dark matter co¬ 
annihilations was proposed by Donato et ah, in 
2000 (H The antideuteron flux at the top of the 
atmosphere due to dark matter co-annihilation and 
decay (called primary flux) can be estimated based on 
dark matter density profiles in the galaxy, dark matter 
co-annihilation and decay channels, hadronization and 
coalescence models, propagation models and solar 
modulation models. The primary antideuteron fluxes 
for different dark matter models are shown in Figure 
The solid blue and purple lines are neutralino and 
gravitino LSPs with 30 GeV 0 and m 2,12 ~ 50 
GeVil respectively, while the dashed green line is for 
LZP with ~ 40 GeV 10. The primary antideuteron 
flux has a relatively flat peak at low energy, E ~ 
0.2 GeV/^. The antideuteron flux due to cosmic-ray 
interactions with the interstellar medium (called sec¬ 
ondary/tertiary flux) is also shown as the solid red line 
in Figure [2l[I2[l3l[Hl. Unlike primary antideuterons, 
collision kinematics suppress the formation of low- 
energy secondary antideuterons. Furthermore, the 
interaction rate is drastically decreased since the flux 


of the cosmic-ray protons follows roughly a power 
law, Fp ~ . Therefore, the primary antideuteron 

flux can be about two orders of magnitude larger than 
the secondary/tertiary antideuteron flux at low energy, 
which can provide clear dark matter signatures. 

The Alpha Magnetic Spectrometer (AMS-02) on 
the International Space Station, launched in 2011, is 
the only currently-operating antideuteron search exper¬ 
iment. AMS-02 probes a higher energy region than 
GAPS, using a different detection technique with a 
completely different background 1201 . GAPS comple¬ 
ments AMS-02 on the antideuteron search, which is 
crucial for rare event searches, while uniquely seeking 
for lower energy and lower background antideuterons. 
GAPS also complements other existing and planned in¬ 
direct searches as well as underground direct dark mat¬ 
ter searches and collider experiments, while distinc¬ 
tively probing a wide variety of dark matter annihilation 
and decay models. The expected GAPS antideuteron 
sensitivity (~ 99% CL) as well as the AMS-02 sensi¬ 
tivity ifTSl [161 dll [m and the current upper limit on the 
antideuteron flux obtained by BESS HU are also shown 
in Figure[^ The details will be discussed in Section 

2. GAPS Project 

The GAPS experiment was proposed in 2002 lITOl . 
and the detection concept discussed below was vali¬ 
dated in an accelerator beam test at KEK, Japan in 2005 
on and the prototype flight (pGAPS) at JAXA/ISAS 
balloon facility in 2012 (15] [T^ [TTl. A GAPS science 
mission is proposed to fly from Antarctica in the austral 
summer of 2019-2020. 

2.7. Detection Concept 

The GAPS detection method involves capturing an¬ 
tiparticles into a target material with subsequent forma¬ 
tion of exotic atoms. A time-of-flight (TOE) system 
measures the velocity (energy) and direction of an in¬ 
coming antiparticle. The antiparticle slows down by the 
dE/dX energy loss and stops in the target material, form¬ 
ing an exotic atom in its excited state. The exotic atom 
de-excites with the emission of Auger electrons as well 
as atomic X-rays ED With known atomic number of 
the target, the Bohr formula for the atomic X-ray en¬ 
ergy uniquely determines the mass of the captured par¬ 
ticle Go). Ultimately, the antiparticle is captured by 
the nucleus in the atom, where it is annihilated with 
the emission of annihilation products, such as pions and 
protons. Since the mean number of pions and protons 
produced by the nuclear annihilation is approximately 
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Figure 2: GAPS detection method: An antiparticle slows down and 
stops in the Si(Li) target, forming an exotic atom. The atomic X- 
rays will be emitted as it de-excites followed by the pion (and proton) 
emission in the nuclear annihilation. 


proportional to the number of antinucleons, the pion and 
proton multiplicities provide an additional discriminant 
to identify the incoming antideutron from other parti¬ 
cles, such as antiprotons. This process is illustrated in 
Figure 



Figure 3: GAPS Antideuteron-antiproton identification technique: (1) 
unique atomic X-rays, (2) pion and proton multiplicity, (3) stopping 
range. 

Antiprotons are the main background in the GAPS 
antideuteron measurement since they can also form ex¬ 
otic atoms and simultaneously emit atomic X-rays and 
annihilation products. Antideuterons, however, are dis¬ 
tinguishable from antiprotons with the atomic X-rays 
and annihilation products as discussed above. Further¬ 
more, the stopping range of antideuterons, which is 
about twice as large as that of antiprotons with the same 


(3, provides an excellent discrimination to distinguish 
antideuterons from antiprotons. Figure [^illustrates the 
GAPS antideuteron-antiproton identification technique 
with atomic X-rays, pion and proton multiplicities, and 
stopping range. The details will be discussed below. 

2.2. Instrumental Design 

The GAPS balloon flight will have several unique 
features. It will be the first balloon flight with a very 
large, pixellated lithium-drifted silicon (Si(Li)) detector 
surrounded by a very large TOF system without a pres¬ 
sure vessel. The Si(Li) detector will be used both as 
a target to form the exotic atom as well as a detector 
for atomic X-rays and annihilation products. There will 
be 10 layers of detectors with each layer composed of 
4 inch diameter, 2.5 mm thick Si(Li) detectors. Each 
Si(Li) detector will be segmented into 4 strips and adja¬ 
cent tracking layers separated by 20 cm will have their 
strips positioned orthogonally. The inner TOF system 
(1.6 m X 1.6 m X 2 m) composed of highly segmented 
plastic scintillators will completely surround the Si(Li) 
layers. The top half of the inner TOF system and Si(Li) 
layers will be covered with plastic scintillators (outer 
TOF system, 3.6 m x 3.6 m x 2 m) and the inner and 
outer TOF planes will be separated by 1 m (see Figure 
1^. This geometry is a natural consequence of the cur¬ 
rent GAPS detection concept, which allows us to track 
and count the number of particles produced in the nu¬ 
clear annihilation, while separately identifying atomic 
X-rays from particle tracks. It also permits direct mea¬ 
surement of particle stopping range in the multi-layer 
geometry. 



Figure 4: GAPS detector design: 10 layers of Si(Li) detectors are 
surrounded by the inner and outer TOF plastic scintillators. Note that 
inner TOF paddles are not shown in the figure in order to show the 
detector layers. Each Si(Li) detector is 4 inch diameter, 2.5 mm thick. 
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3. GAPS Antideuteron Sensitivity 

As discussed above, antideuterons can be identified 
from other cosmic-ray particles by using the GAPS de¬ 
tection technique with atomic X-rays and pion/proton 
multiplicities, as well as the different stopping range 
of incoming particles. Low-energy cosmic-ray protons, 
even though their flux is several orders of magnitude 
larger than the antideuteron flux, can be clearly dis¬ 
tinguishable from antideuterons since they cannot emit 
atomic X-rays nor annihilation products. Furthermore, 
the coincidence of multiple cosmic-ray protons coming 
into the same channel within a short time window (~ 
100 ns) is negligible. Therefore, we will focus on an¬ 
tiprotons as the main background in the GAPS exper¬ 
iment as they can also form exotic atoms and produce 
atomic X-rays and charged particles. In the following 
section, we will discuss how antideuterons can be iden¬ 
tified from antiprotons and also show the simple, but 
comprehensive approach to the GAPS antideuteron sen¬ 
sitivity estimation by using Monte Carlo simulations, 
atomic cascade model and Intra-Nuclear Cascade (INC) 
model. The GAPS antideuteron sensitivity was esti¬ 
mated with -99% confidence level (CL) to identify pri¬ 
mary antideuterons from other cosmic-ray particles. 

3.1. Grasp 

Antiparticles coming from the top of the atmosphere 
are degraded and also annihilated in-flight in the atmo¬ 
sphere before reaching the GAPS instrument. Some of 
them are able to reach and stop inside the GAPS instru¬ 
ment. A Geant4 Monte Carlo simulation was conducted 
to evaluate how many antiparticles can reach the GAPS 
instrument during the balloon flight. As seen in Figure 

antiprotons and antideuterons were generated above 
the compressed atmosphere (a fiat plane, column den¬ 
sity - 4g/cm^) to evaluate the energy loss and in-flight 
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Figure 5: The Geant4 setup for grasps estimation. Antiprotons and 
antideuterons were generated above the compressed atmosphere (a flat 
plane, column density ~ 4g/cm^) to evaluate the energy loss and in¬ 
flight annihilation in the atmosphere. 


annihilation in the atmosphere. The grasps [rn^sr] of 
the GAPS instrument at - 40 km flight altitude for an¬ 
tiprotons and antideuterons that stop in the Si(Li) layer 
Ystop annihilate in-flight in the instrument (not in 
the atmosphere) F^"^ were estimated with Geant4.10.01. 
Note that F^^^^ is defined as a product of the effective 
area, solid angle and stopping efficiency, while F^"^ is 
defined as a product of the effective area, solid angle 
and efficiency of the in-flight annihilation. 



Figure 6: The grasps for antiprotons (red lines) and antideuterons 
(black lines). The solid lines are grasps for stopped events, while the 
dashed lines are for in-flight annihilation events 


The grasps for antiprotons (black lines) and an¬ 
tideuterons (red lines) in terms of the energy at the top 
of the atmosphere are shown in Figurej^ The solid lines 
are grasps for stopped events, while the dashed lines are 
grasps for in-flight annihilation events. Note that sepa¬ 
rate calculations show other antideuteron loss processes, 
such as Coulomb breakup or the Oppenheimer-Phillips 
effect, to be small (< 10%) 1221 . 


3.2. Atomic X-ray 

As discussed above, the stopped antiparticles can 
form exotic atoms and emit atomic X-rays during the 
de-excitation. The energies of the atomic X-rays are 
uniquely determined by the components of the exotic 
atom as seen below. 


M* 




1 1 


(A) 


Here, z and Z are the charge of the antiparticle and tar¬ 
get atom, M* and m* are the reduced masses of an an¬ 
tiparticle in the exotic atom and an electron in the tar¬ 
get atom, Rh is the Rydberg constant and Hi and rif are 
the initial and final principal quantum number. The an- 
tideuteronic exotic atom formed with the Si target can 
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emit 30 keV, 44 keV and 67 keV atomic X-rays, while 
23 keV, 35 keV and 58 keV X-rays can be emitted for 
antiprotonic exotic atoms formed with the Si target. The 
X-ray yield, defined as the probability to emit atomic 
X-rays per exotic atom, for each atomic X-ray was esti¬ 
mated as ~ 80% with the simple cascade model as dis¬ 
cussed in ED- Note that in addition to the antiparticles, 
atomic X-rays were manually generated in the simula¬ 
tion since the cascade model of the exotic atoms has 
not been implemented yet in the current Geant4 physics 
package. Figure [7] shows the simulation results for en¬ 
ergy spectra in the Si(Li) detector, taking into account 
all the physics processes and interactions of particles 
and radiation passing through matter. The black and red 
lines are for antiproton and antideuteron events, respec¬ 
tively. The dashed lines show energy spectra without 
detector response, while the solid lines taking into ac¬ 
count the FWHM energy resolution (~ 3 keV) of the 
Si(Li) detector. The continuous spectra are due to the 
EM shower developed by the annihilation products. The 
antideuteron selection cuts and efficiencies to atomic X- 
ray detection are shown in Table The efficiencies to 
the detection of at least one antideuteronic atomic X-ray 
with these selection cuts were estimated as ~ 0.02 
for antiprotons and ~ 0.12 for antideuteron events, 
as below. 



Figure 7: Geant4 simulation results for energy spectra in the Si(Li) de¬ 
tector. The black and red lines are antiproton and antideuteron events, 
respectively. The dashed lines show energy spectra without detector 
response, while the solid lines taking into account the FWHM energy 
resolution (~ 3 keV) of the Si(Li) detector. 


^ X >\ - 1 - (1 - 63 o )(1 - 644)(1 - 657) 

Here, 630 , 644 and 657 are the efficiencies to an¬ 
tideuteronic atomic X-rays with the selection cuts. Note 


that atomic X-rays for antideuterons stopped in the de¬ 
tector frame and other materials are also distinguishable 
from those for antiprotons, since their atomic X-rays 
are different with each other, as seen in Eq Q. The 
stopped position of the incoming antiparticle can be pre¬ 
cisely determined by tracking and timing the annihila¬ 
tion products (charged pions) backwards in the detector 
layers and the TOE paddles. 

Table 1: The select cuts and efficiencies to the antideuteronic atomic 
X-rays for antiproton and antideuteron events 


Energy 

selection cuts 



30keV 

29.5 <E< 30.5 keV 

9.5 xlO-^ 

3.8 xlO-2 

44keV 

43 < E < 44.5 keV 

8 . 0 x 10 -^ 

5.4x10-2 

67keV 

66 < E < 68 keV 

6.7x10-^ 

3.0x10-2 


3.3. Nuclear Annihilation Products 

3.3.1. Intra-Nucear Cascade (INC) Model 

The interaction of antiprotons with nucleons (pro¬ 
tons and neutrons) has been studied since the dis¬ 
covery of the antiproton in 1955 Il23]| . The Intra- 
Nuclear Cascade (INC) model has been developed to 
predict the particle multiplicity in the antiproton anni¬ 
hilation on nuclei GSI^. The INC model is com¬ 
posed of four stages: ( 1 ) primordial pion production, 
( 2 ) direct emission (fast ejectiles) from the primordial 
pion-nucleon interaction, (3) pre-equilibrium emission 
(multi-fragmentation) from excited nucleus, and (4) nu¬ 
clear evaporation (slow ejectiles). The antiproton is 
first assumed to annihilate at the surface of a single nu¬ 
cleon of the nucleus and to produce the primordial pi¬ 
ons. Some of these pions may escape from the nucleus, 
but the others cascade through the nucleons of the nu¬ 
cleus, knocking out fast nucleons as they go through. 
Due to this interaction, the nucleons in the nucleus are 
excited and their energy density is far from the thermal 
equilibrium state. Therefore, the nucleus can break into 
pieces with the emission of fragmented particles. The 
characteristic time of this process is ~ 10“^^ s. Then 
the density distribution becomes more thermal and the 
remaining excitation energy will be removed by nuclear 
evaporation, which emits slow ejectiles. Figureshows 
an overview of the INC model processes. 

The INC model can also be applied to the an¬ 
tideuteron annihilation to estimate pion and proton mul¬ 
tiplicities, resulting in two models that have been built 
based on how the two antinucleons in the antideuteron 
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Figure 8: Schematic of INC model. (1) The primordial pious are 
produced in the nuclear annihilation. (2) Some of the pions hit the nu¬ 
cleons in the nucleus with direct emission (fast ejectiles), followed by 
(3) the pre-equilibrium emission (fragmentation) and (4) the nuclear 
evaporation (slow ejectiles). 



Figure 9: Schematic of the interaction probability for primordial pions 
with the nucleons in the atom 


interact with nuclei 1^ . The first model (model A) as¬ 
sumes that the two antinucleons interact with the nu¬ 
cleons simultaneously, while the second model (model 
B) assumes that the antinucleons interact with nucleons 
separately. We will discuss how we can use the INC 
model to distinguish antideuterons from antiprotons in 
the following sections. The antideuteron selection cut 
and the efficiency were numerically evaluated. 


3.3.2. Pion Multiplicity 

The number of primordial charged and neutral pi¬ 
ons in the INC model is estimated based on antiproton- 
nucleon annihilation. This has been well studied and 
much experimental data are available. The pion multi¬ 
plicity and the standard deviation (cr) for the antiproton- 
nucleon annihilation is estimated as follows 1 ^1^ : 


(M^,o> = 2.65-f 3.65 In 




= 0.174 (Vs) 


0.40 


(B) 

(C) 


Here, is the average number of primordial 

charged and neutral pions and ^/s is the center of mass 
energy in GeV. Eq ^ and Eq ^ agree well with the 
experimental data ||24l[27l[28l. The average number of 
primordial charged and neutral pions and the av¬ 

erage number of primordial charged pions ) for the 
antiproton-nucleon annihilation at rest are 5.1 and 3.1, 
respectively. 

Most of the primordial pions can escape from the nu¬ 
cleus, but some of them interact and cascade through 
the nucleons, producing emission of fast nucleons. The 


interaction probability can be geometrically calculated 
for each atom A, and the average number of the (fi¬ 
nal) charged pion multiplicity {M „±) can be simply esti¬ 
mated as follows 12 ^1^ : 


P(A) 

Q(A) 


<M^,). P(A) 
fl(A) 


1 - 


471 


27r(l - cosO) 
( 


271 


L ( rpA^^ f 

\ \roAi/3+(5/ 


Here, P(A) is the pion-nucleus interaction probability, 
Q(A) is the solid angle that the pion can hit the nucleons 
in the atom, and rq and S are the radius parameters, 1.2 
fm and 1.6 fm, respectively |[30l. The INC model agrees 
well with the antiprotonic experimental data 1^ . Note 
that the pion multiplicity reduces by only 0.3 as the nu¬ 
clei changes from A = 10 to A = 200. 

The primordial charged and neutral pion multiplicity 
for the antideuteron annihilation at rest on nuclei was 
estimated with Eq and Eq ^ by simply changing 
= 4 GeV in model A and ^[s = 2 GeV for each 
antinucleon in model B. The (final) charged pion mul¬ 
tiplicity can be estimated with the pion-nucleus inter¬ 
action probability as discussed above. The pion multi¬ 
plicity and the efficiency to selection cuts for antipro¬ 
ton (6^) and antideuteron (6^) annihilations at rest on 
the Si target are shown in Table Here, we assume 
{Mj^+) = = {Mj^q) and consider the average pion 


6 










multiplicities between model A and B. The estimated 
pion multiplicity and the efficiency to selection cuts for 
in-flight annihilation events (annihilation with the ki¬ 
netic energy T ~ 0.1 GeV/n) were also shown in Table 

m 

Table 2: The pion multiplicity and the efficiency to selection cuts for 
antiproton and antideuteron annihilations at rest and in-flight (annihi¬ 
lation with the kinetic energy T ~ 0.1 GeV/n) on the Si target. 


however, estimated the nucleon (proton -l- neutron) mul¬ 
tiplicity with E >60 MeV for the antideuteron annihila¬ 
tion at rest on the Mo (Molybdenum) target. The proton 
multiplicity for the Si target was simply estimated by 
scaling the multiplicity of the Mo target to the Si tar¬ 
get by taking into account the pion-nucleus interaction 
probability as discussed above. We also assumed the re¬ 
lationship between the proton and neutron multiplicity 
as follows 1291 



annihilation at rest 

in-flight 

{M„.) 


4 

4 

4 

>3 

0.53 

0.93 

0.57 

0.93 

> 4 

0.25 

0.80 

0.58 

0.82 

>5 

0.070 

0.62 

0.092 

0.64 

>6 

0.093 

0.41 

0.019 

0.44 

>7 

0.00067 

0.23 

0.0024 

0.26 

>8 

0.000026 

0.11 

0.00016 

0.13 


3.3.3. Proton Multiplicity 

The INC model also predicts the proton and neutron 
production in the following processes: (1) direct emis¬ 
sion from the interaction between the primordial pions 
and the nucleus, (2) pre-equilibrium emission (multi 
fragmentation) from excited nucleons, and (3) nuclear 
evaporation. The direct emission process dominates for 
protons with E > 30 MeV and the energy spectrum for 
the proton multiplicity is estimated with the Maxwell- 
Boltzmann distribution as given below, based on fitting 
of the experimental data for antiproton annihilation at 
rest on nuclei 12 ^ : 


{Mn)l{Mp) ^ 

NjZ 

Here, N is the number of neutrons in the nucleus, and Z 
is the atomic number. Considering the above, the pro¬ 
ton multiplicity with E >60 MeV for the antideuteron 
annihilation on the Si atom was estimated as 2.35 for 
model A and 1.76 for model B. Table shows the pro¬ 
ton multiplicity with E > 60 MeV and the efficiency 
to selection cuts for antiproton and antideuteron (aver¬ 
age between model A and B) annihilations at rest on the 
Si atom. Here, the proton multiplicities for antiprotons 
and antideuterons were simply estimated as the Poisson 
distribution. The proton multiplicity and the efficiency 
to selection cuts for in-flight annihilation events (anni¬ 
hilation with the kinetic energy T ~ 0.1 GeV/n) were 
also estimated simply by scaling the proton multiplicity 
based on the primordial pion multiplicity as seen in Eq 

Table 3: The proton multiplicity with E >60 MeV and the efficiency 
to selection cuts for antiproton and antideuteron annihilations at rest 
and in-flight (annihilation with the kinetic energy T ~ 0.1 GeV/n) on 
the Si atom. 


dM 



Here, (M) is the average number of protons produced in 
the antiproton annihilation at rest on Si, 0.86 + 0.05, E 
is the energy of the proton, and T is the parameter fitted 
to the data (~ 40 MeV). If the energy of the proton is 
too small, it will stop quickly and could not be tracked 
by the detector layers. Therefore, we set the lower cut 
of the proton energy as 60 MeV to guarantee passage 
through three or more Si(Li) layers. The proton multi¬ 
plicity with E >60 MeV <M^^) is 0.37. Note that these 
protons are distinguishable from pions generated in the 
nuclear annihilation since, unlike pions, they are not rel¬ 
ativistic and deposit much more energies in the detector 
layer. 

Unfortunately, there is no data available for the an¬ 
tideuteron annihilation at rest on Si. Cugnon et al. ll^ . 



annihilation at rest 

in-flight 

{Mf) 

4 

P 

4 

P 

> 1 

0.31 

0.87 

0.32 

0.88 

>2 

0.054 

0.60 

0.057 

0.62 

>3 

0.0064 

0.34 

0.0070 

0.36 

>4 

0.00058 

0.16 

0.00066 

0.0.16 


3.4. Depth Sensing and dEjdX Energy Deposit 

The incoming particle can be tracked in the TOE 
paddles and Si(Li) detectors, as discussed above. The 
number of layers that the incoming particle can pene¬ 
trate before stopping in the detector, provides the stop¬ 
ping range (called depth sensing) of the incoming parti¬ 
cle. The stopping range for the antideuteron is approxi¬ 
mately twice as large as that for the antiproton with the 
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same velocity, as seen in the scaling law of the stopping 
range below uni. 






mp z 


2^P I 


m3 




Here, Rj(Rp) is the stopping range of antideuterons (an¬ 
tiprotons), (mp) is the mass of antideuterons (an¬ 
tiprotons), Zd {zp) is the charge of antideuterons (an¬ 
tiprotons) and Tj (Tp) is the kinetic energy of an¬ 
tideuterons (antiprotons). Therefore, antideuterons can 
be distinguished from antiprotons by using the depth 
sensing and the TOP timing. Note that the stopped posi¬ 
tion of the incoming antiparticle can be precisely deter¬ 
mined by tracking and timing the annihilation products 
(charged pions) backwards in the detector layers and the 
TOP paddles, as discussed above. 

A Geant4 simulation was conducted to evaluate the 
depth sensing for antideuteron selection cuts. Protons 
and deuterons were used to simply evaluate the depth 
sensing for antiprotons and antideuterons. The simula¬ 
tions were done for two different incoming angles, 0 deg 
and 45 deg, taking into account the timing and angular 
resolution of the TOP system, 0.5 ns and ~ 5 deg. The 
tracking capability allows us to identify what material 
the incoming particle went through before stopping. 


Antideuterons stopped in the same layer as antiprotons 
show a larger TOP timing, as expected. The efficiency 
to the selection cut, tpoF > 13.3 ns (17.6 ns) for antipar¬ 
ticles coming with 0 deg (45 deg) becomes 5.3 x 10“"^ 
(6.1 X 10“"^) for antiprotons, while 0.83 (0.78) for an¬ 
tideuterons. 



TOP [ns] 


Figure 11: The TOF timing for antiprotons (black lines) and an¬ 
tideuterons (red lines) with the incoming angle 45 deg and stopped 
in layer 3. The dashed lines show the TOF timing without detector re¬ 
sponse, while the solid lines taking into account the timing resolution 
(0.5 ns) of the TOF system 



Table [ 4 ] illustrates the antideuteron selection cuts and 
corresponding efficiencies for particles stopping in the 
different Si(Li) layers. The depth sensing provides 
excellent antideuteron identification from antiprotons, 
4epth = 6-0 X 10-4 and 4^,^ = 0.68 on average. Note 
that dP/dX energy deposit in the Si(Li) detector can be 
also used to distinguish antideuterons from antiprotons 
since antideuterons with the same TOP as antiprotons 
can have roughly twice as much kinetic energy as an¬ 
tiprotons, and deposit more energy in the Si(Li) detec¬ 
tor. This will be evaluated with more detailed simula¬ 
tions and analyses, such as an event-by-event analysis 
and likelihood analysis, in the future. 


Figure 10: The TOF timing for antiprotons (black lines) and an¬ 
tideuterons (red lines) with the incoming angle 0 deg and stopped in 
layer 3. The dashed lines show the TOF timing without detector re¬ 
sponse, while the solid lines taking into account the timing resolution 
(0.5 ns) of the TOF system 

Pigure p^ ( pT] ) shows the TOP timing for antiprotons 
(black lines) and antideuterons (red lines) with the in¬ 
coming angle of 0 deg (45 deg) and stopped in layer 
3. The dashed lines show the TOP timing without de¬ 
tector response, while the solid lines taking into ac¬ 
count the timing resolution (0.5 ns) of the TOP system. 


3.5. Confidence Level and Sensitivity Calcuration 

The major background events in the GAPS an¬ 
tideuteron search are antiprotons that pass all the an¬ 
tideuteron selection cuts, as well as the secondary an¬ 
tideuterons produced by cosmic-ray interactions. The 
antiproton and the secondary antideuteron background 
events can be estimated as follows. 

= f Fp(E)rp(E)T€^eTdE ■ ef 
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Table 4: The antideuteron select cuts and the efficiency to the depth 
sensing for antiproton and antideuteron events. The incoming angles 
for antiparticles are 0 + 5 deg and 45 + 5 deg. 


exotic atoms. The total antideuteron sensitivity with ~ 
99% CL, was estimated as below. 

a 



0 deg 

45 deg 

layer 

cuts 

depth 

depth 

cuts 

eP 

depth 

depth 

1 

> 16.7 ns 

5.9x10-4 

0.66 

> 24.2 ns 

7.2 xlO-4 

0.82 

2 

> 14.5 ns 

5.2x10-4 

0.81 

> 19.4 ns 

6.9x10-4 

0.61 

3 

> 13.3 ns 

5.3x10-4 

0.83 

> 17.6 ns 

6.1 xlO-4 

0.78 

4 

> 12.5 ns 

6.1 xlO-4 

0.81 

> 16.7 ns 

4.9 xlO-4 

0.75 

5 

> 12.0 ns 

5.4x10-4 

0.74 

> 15.6 ns 

4.6 xlO-4 

0.79 

6 

> 11.5 ns 

6.7 xlO-4 

0.72 

> 15.0 ns 

7.8 xlO-4 

0.83 

7 

> 11.1 ns 

7.7x10-4 

0.70 

> 14.6 ns 

6.3 xlO-4 

0.77 

8 

> 10.8 ns 

7.5x10-4 

0.66 

> 14.2 ns 

7.4x10-4 

0.73 

9 

> 10.7 ns 

4.7x10-4 

0.51 

> 13.8 ns 

7.5 xlO-4 

0.72 

10 

> 10.4 ns 

5.2x10-4 

0.55 

> 13.7 ns 

1.3 xlO-4 

0.55 

ave 

- 

6.0x10-4 

0.70 

- 

6.0 xlO-4 

0.66 


^bkg ^ J Fi^^(E)rj(E)T€^€TdE ■ ef 

Here, Fp and are the antiproton and secondary an¬ 
tideuteron fluxes at the top of atmosphere, Fp and Fj 
are the grasps for antiprotons and antideuterons, 6 ^ is 
the geomagnetic cutoff parameter for antideuterons, T 
is the flight time, ej is the ratio of the observation time 
to the flight time (~ 1 . 0 ), ef is the efficiency to each 
cut and i is the cut type (atomic X-ray, pion/proton mul¬ 
tiplicity and depth sensing). Note that GAPS is pro¬ 
posed to fly from Antarctica, which is the best loca¬ 
tion with the lowest geomagnetic cutoff. Eg ~ 0.45 for 
antiprotons and ef ~ 0.75 for antideuterons at ~ 0.2 
GeV/^ (evaluated with the Geant4-based PLANETO- 
COSMICS framework |[3^ [3^ ). The sensitivity to n 
antideuteron detections with the selection cuts, 5 5 , and 

a 

the corresponding confidence level (CL) to identify an¬ 
tideuterons from background events can be estimated as 
follows. 


an 


CL 


1 


\^F^{E)TeleTdE^Y\ief 
\-P{N = + 


P{N,X) = 


Nl 


Here, P(N,A) is the Poisson distribution with the ob¬ 
servation value N and the mean A. The selection cuts 
for stopped events and in-flight annihilation events are 
different since in-flight annihilation events do not form 


total 

d 


V d,stop 


+ 


d,inj / 


(D) 


Here, the antideuteron selection cuts for stopped events 
are ’’depth sensing -\-7i>5 -f (X>i or p >2 or 7 r> 2 )” for three 
LDB flights (LDB, 105 days), while they are ”7r>6 + 
(p >2 or 7 r> 2 )” for in-flight annihilation events. Two or 
more antideuteron detection is required for in-flight an¬ 
nihilation events to obtain -99% CL. The correspond¬ 
ing antideuteron sensitivity and the number of back¬ 
ground events are 2.0 xl 0 “^ [m“^s“^sr“^(GeV/^)“^] and 
~ 0.01 events, respectively. Note that the number of 
background events increases to 0.014 if atmospherically 
produced antiprotons and antideuterons are taken into 
account ca. Atmospheric antiproton and antideuteron 
production and transport through the atmosphere and 
geomagnetic held will be studied more in the future. 

The GAPS antideuteron sensitivity as well as an¬ 
tideuteron fluxes at the top of the atmosphere for dif¬ 
ferent dark matter models are shown in Figure Here, 
GAPS is sensitive to neutralino and gravitino LSPs in 
SUSY theories, as well as LZP in extra dimension the¬ 
ories, discussed in Section |1.2| Moreover, the GAPS 
antideuteron sensitivity is two orders of magnitude bet¬ 
ter than the current upper limit on the antideuteron flux 
obtained by the BESS experiment |[T9]| . 

The AMS-02 antideuteron sensitivity (five years of 
observation time) was estimated for two different en¬ 
ergy regions where two different detectors, the time of 
flight (0.2 GeV/^ < E < 0.8 GeV/^) and the ring imag¬ 
ing Cherenkov counter (2.2 GeV/^ < E < 4.2 GeV/^), 
are used for the velocity measurement (1811201. Here, 
the geomagnetic cutoff at the International Space Sta¬ 
tion was taken into account for the sensitivity calcula¬ 
tion. The antideuteron sensitivity in 2.2 GeV/^ < E < 
3.6 GoVIn is not as sensitive as the one in 3.6 GoVIn < 
E < 4.2 GoVIn, due to the difference of the antideuteron 
identification capability. These two regions were not 
combined using Eq (|^. This was done to maintain con¬ 
sistency with the AMS-02 sensitivity calculations of 
dsmainiiii. All the sensitivities were estimated for 
-99% CL to identify primary antideuterons from other 
cosmic-ray particles. Note that the AMS-02 sensitivity 
calculation is optimistic since it is based on the super¬ 
conducting magnet, rather than the permanent magnet 
used in the actual flight. 
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4. Conclusion 

Cosmic-ray antideuterons can provide a clean dark 
matter signature since the primary antideuteron flux due 
to the dark matter co-annihilation and decay can be 
a few orders of magnitude larger than the secondary 
flux at low energy. GAPS can uniquely probe dark 
matter through low-energy antideuterons with distinc¬ 
tive detection methods, while complementing exist¬ 
ing and planned other indirect dark matter searches as 
well as underground direct dark matter searches ex¬ 
periments and collider experiments. The simultane¬ 
ous detection of atomic X-rays and charged particles 
from the decay of exotic atom as well as the timing 
and depth sensing of the incoming particle provides ex¬ 
cellent particle identiflcation capabilities to select an¬ 
tideuterons from other cosmic-ray particles. The GAPS 
antideuteron sensitivity was estimated to be 2.0 xl0“^ 
[m“^s“^sr“^(GeV/^)“^] for three LDB flights (105 days, 
~ 99% CL), which is more than two orders of magnitude 
better than the current upper limit on the flux obtained 
by BESS, and also sensitive to various dark matter mod¬ 
els. 
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